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Abstract We have computed a grid of chemical evolution 
models for a large set of spiral and irregular theoretical galax- 
ies of different total mass. In our models, the gas phase has 
two components, the diffuse and the molecular one (H2). It 
is possible, therefore, to follow the time (or redshift) evolu- 
tion of the expected density of the H2 phase. We will show 
the predictions of this gas density at hight redshift, which 
might be detected with ALMA, in this type of galaxies. 

Keywords Galaxies:Spirals • Galaxies: Irregulars • 
Interstellar medium • Molecular gas 

1 Introduction 

The molecular interstellar medium plays a critical role in 
the evolution of galaxies since it provides the material from 
which stars form. There are observations [9) of rotational 
transitions of CO within 1 .0 < z < 6.5, mostly in quasars and 
Early Universe Molecular Emission Line Galaxies (here- 
inafter QSO's and EMGs, respectively) which demonstrate 
that molecular clouds exist already at the epoch of the galaxy 
formation. However, the current sample is scarce for z > 3, 
where the study of the star and galaxy formation is spe- 
cially interesting. ALMA is the instrument that may solve 
this problem. 
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In our chemical evolution models the gas phase has two 
components, the diffuse and the molecular phases of gas, 
which we treat separately. This way, it is possible to follow 
the evolution of the expected density of H2 with redshift. We 
will show here our predictions about the H2 density at high 
redshift, and we will analyze if they might be detectable with 
ALMA. 



2 A grid of multiphase chemical evolution models 

We have computed a wide grid of chemical evolution mod- 
els for a large set of theoretical spiral and irregular galaxies 
of different total mass [5] by using the multiphase chemical 
evolution model from [3 4 6]. We have simulated theoreti- 
cal models with 44 different radial distributions of total mass 
and 10 star-forming efficiencies, with values in the range [1- 
0], which we denoted by N from 1 to 10, for each one of 
them. The 440 resulting models provide the time evolution 
for specific radial regions within each galaxy modeled, with 
the corresponding mass in each gaseous or stellar phase, as 
well as their elemental abundances. This type of model was 
already applied to several individual spiral galaxies in Q, 
where it was shown that the model is adequately calibrated 
and that it may be used for galaxies of different morpholog- 
ical types and different total masses. 

The code assumes that a total mass, M„ a \, initially as dif- 
fuse gas, is located within a protogalaxy or halo. The gas 
falls onto the equatorial plane forming out the spiral disk 
at a rate defined by a collapse time-scale which depends on 
M ga i. In each zone, halo or disk, several radial regions are 
defined. We follow separately the evolution in each region. 

The assumed star formation in the halo is a power law of 
the diffuse gas. In the disk, however, we assume that the star 
formation occurs in two steps: first, molecular clouds (H2) 
form from the diffuse gas (HI) also following a power law. 
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Fig. 1 Radial distributions of molecular gas surface density at z = 0. Each panel represents the results for a different rotation velocity galaxy, 
which means a different galaxy total mass. Curves are for different star formation efficiencies in the range 1 (the highest efficiency, N = 1) to 
(the lowest efficiency, N = 10) as labelled in panel a). 



Then, stars form by two means: 1) cloud-cloud collisions, 
and 2) by the interaction of the massive stars with the sur- 
rounding molecular clouds: This way, in our models the gas 
has two components: HI and H2. 

Since we have computed separately the diffuse gas evo- 
lution and the molecular clouds evolution, we obtain the ra- 
dial distributions of the molecular gas mass, which is repre- 
sented in Fig.Q] There we show the radial distributions of the 
molecular gas density at the present-time for 4 values of the 
rotation velocity, or total mass of the Galaxy. Thus, panel a) 
shows the radial distributions produced in a low mass galaxy 
while d) represents a massive galaxy. Panel c) is the modeled 
galaxy most similar to our Milky Way Galaxy and panel d) 
is an intermediate mass galaxy. In each one of these panels 
models for different star formation efficiencies, named by N, 
are represented. The lines in color correspond to the lowest 
star-forming efficiencies (high N values), which maintain a 
higher density than in cases where stars are formed more 
efficiently (black lines). The distributions are close to ex- 
ponential in the external regions but are flatter at the inner 
regions, mostly when the N is low, as is in fact observed in 
Nature 0. 



The evolution of the total mass of H2 in galaxies is rep- 
resented in Fig|2] The magenta open stars are estimates from 
1 1 1 for nearby QSOs. The predicted present time molecular 
gas masses are similar to those ones of local objects with 
a maximum for z — in agreement with the observational 
limit log(MH2) = 10.4, shown by a green solid line, ob- 
tained for those authors. The red full dots are the estimates 
from [9 1 for EMGs. The observed molecular gas masses are 
decreasing for decreasing z in agreement with the most of 
our model results. Only the models with low star-forming 
efficiencies (colored lines) are still producing H2. 

On the other hand, following that work |9|, and taking 
into account the estimates for the star formation rate and 
molecular gas total mass, these objects might be bulges or 
starbursts in centrally concentrated disks, but not elliptical 
galaxies as it was initially assumed. Our models are also in 
agreement with this statement: the total mass of H2 for the 
most massive spiral galaxies (with Vrot in the range 100-300 
km.s ) and high star formation efficiencies (N < 5) are al- 
ways below of the EMGs data region since in those cases the 
molecular gas is very quickly consumed. Only the models 
for massive galaxies -panels c) and d) with the lowest effi- 
ciencies show values closer to the observations of EMG's. 
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Fig. 2 The evolution of the total mass of molecular gas for some of our models. Each panel represents the results for a given galaxy total mass (or 
Vrot), using, for each one, 10 different SFR efficiencies. Lines mean the same as in Fig.l The black lines correspond to the highest efficiencies. 



At this high redshift, these particular models show star for- 
mation only at the central regions and therefore they would 
appear as concentrated disks, with a large quantity of H2 . 

3 The CO detection limits 

Observation of emission from CO rotational transitions is 
the usual means of detecting interstellar molecular clouds. 
CO is a very stable molecule and the most abundant after 
H2. Following [2 1 levels J = 4 — 5 are the best ones to detect 
molecular gas at high redshifts. 

In order to verify if our densities of molecular gas may 
be detected with ALMA, we need to transform the limiting 
flux for this instrument to a limit in gas density. The flux 
limit depends on the observed frequency. We take the high- 
est value from the report "Science with ALMA", and we use 
Eq. 1 from 0, to calculate the minimum Leo necessary for 
detection. This luminosity is transformed into gas mass by 
using, as those authors, a constant conversion factor a = 0.8 
(which assumes independence of redshift z). The result is 
shown as a thick line in Fig [3] where we show the expected 
radial distributions of molecular gas density normalized to 
a characteristic radius Rear, taken as the half optical radius, 
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Fig. 3 The expected radial distributions of molecular gas density nor- 
malized to Rear (half optical radius) for three redshifts as labelled. 



for each total mass radial distribution. We show in top panel 
a) the results for z — (present time) while panels b) and c) 
represents the same models for the higher redshifts z = 2.95 
(t ~ 1.2 Gyr) and z — 4.25 (t ~ 0.6 Gyr), respectively. For 
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Fig. 4 Expected evolution of the flux density of CO lines as a function of z for spiral galaxies with four rotation velocities. Line meaning as in 
Fig.l 



z = most modeled galaxies show distributions with values 
higher than the limiting value which may be detected with 
ALMA down to the optical radius, while for high redshifts 
only our inner disks, mainly of low efficiencies (colored) 
models, have a flux density sufficiently high to be detected. 
Therefore, if the red points of Fig. 2 would correspond to 
spiral galaxies, they must represent, effectively, their central 
parts. 

In Fig. |4] we represent the expected evolution of the flux 
density of CO lines as a function of z, for an integrated 
galaxy, computed from the molecular gas density cr(H2) ob- 
tained from our models, transformed to fluxes by using the 
method outlined. The most massive galaxies -panel d- and 
those with high star formation efficiencies -black lines- main- 
tain an almost constant H2 density and flux density up to 
Z ~ 10, with values just over the detection limits. Thus, such 
as we have before claimed, the most efficient galaxies in 
forming stars, which consume very quickly the molecular 
gas, are not the most probable objects to be detected. Those 
with low rotation velocities Vrot < 50 km.s-1, on the other 
hand, show increasing densities for decreasing z, having val- 
ues lower than the limiting flux for z > 0.5. Only the massive 
galaxies Vrot > 200 km.s -1 with N > 5 have fluxes above 
the observational limit even at z ~ 8. Therefore, the most 



plausible objects to be observed with ALMA will be the cen- 
ters of low-intermediate efficiencies massive galaxies. 

4 Conclusions 

- At redshift z > 3 only the central regions of galaxies 
might be observed in CO with ALMA 

- EMGs are not necessarily the largest or the most massive 
stellar disks. 
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